Generalized Schwinger boson mean field theory for two-dimensional SU(4) symmetric 

spin-orbital antiferromagnet 
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A generalized Schwinger boson representation is proposed to describe the two-dimensional SU(4) 
symmetric spin-orbital quantum antiferromagnet. A uniform mean field solution gives rise to an 
antiferromagnetic long range ordered state with a small staggered magnetization m s = 0.07566, 
and the recent quantum Monte Carlo simulations on rather large lattice sizes have provided strong 
evidence to support these mean field results. 

PACS numbers: 75.10.Jm 



There have been much interests in the properties of 
Mott insulators with orbital degeneracy. In order to de- 
scribe low-energy physics of an insulating phase with one- 
electron per site with double orbital degeneracy in the 
limit of the large Hubbard repulsion, an SU(4) symmet- 
ric model was proposed [l| 



H 



<i,j> 



[S i -S i +T i -T i -4(S i -S i )(T i -T J -)] ! 



where Si and Tj denote the spin- 1/2 and orbital- 1/2 
operators at a lattice site i, respectively. The physi- 
cal condition to derive the above effective model is that 
among the possible two-particle states obtained upon vir- 
tual hopping, the inter-orbital singlet is the lowest energy 
state due to the dynamic Jahn- Teller effect. 

Recently, by generalizing the SU(2) Schwinger boson 
representation |2j, we proposed an SU(4) spin-boson rep- 
resentation to denote both the spin-1/2 and orbital-1/2 
operators simultaneously, and the model is simplified to 
a quadratic form of a nearest neighboring SU(4) symmet- 
ric valence bond (VB) pairing operator |3j. By a direct 
Hartree-Fock decomposition and replacing the local con- 
straint by its bulk average, a uniform SU(4) mean field 
(MF) theory gives rise to a spin-orbital liquid state with a 
finite gap in ID, consistent with the exact matrix product 
solution Q and the numerical simulations pj. On a 2D 
square lattice, the spin, orbital, and spin-orbital tensor 
form antiferromagnetic (AF) long range ordering state. 

It is noted that total spin, orbital, and staggered spin- 
orbital operators Ej 5 ". T,j If, V 2V.>7 /< H with 
Q the AF reciprocal vector, generate an SU(4) Lie alge- 
bra and commutes with the Hamiltonian. In order to 
characterize both spin and orbital degrees of freedom 
simultaneously, an SU(4) generalized Schwinger boson 
representation is introduced as the generators Fg(i) — 

b\ Jai.fi with the expressions 



bt A h 



4 u i,3, 



S* = (b\ A b iA - b[ 2 b t , 2 + b\ 3 b h3 - b\ A b iA )/2; 
T+ = < As + b\ 2 b iA , 77 = 6] i3 6 M + 6t 4 6i, a , 
Tt = (<A,i + J>J,Aa - b{ 3 b i}3 - b\ A b iA )/2, (1) 



with a local constraint Eu^I/Am = The m °del 
Hamiltonian is then rewritten as 



H = ~ E B h B u 



(2) 



<i,]> 



where B it j = [&i,A,i + A,4 - &j,3&i,2 - &i,2&i,3]) an 
SU(4) symmetric singlet pairing operator. 

When a VB order parameter is assumed A = 
— (Bi.j) and the local constraint is imposed by a La- 
grangian multiplier A, in terms of a Nambu spinor 



k.2! u k,3' u k,4' b-k,l, &-k,2, &-k,3, &-k,4j 



the MF Hamiltonian matrix is obtained as H m f(k) = 
A + 2ZA 7k f2i, where Q\ — a x ® <J y ® a yi and 7 k 



1 E<5 



k<5 



From the Lagrangian, the bosonic Matsub- 



ara Green function is given by 



G(k,zw r . 



—iw n tt2 — A + 2ZA7 k 51i 
W 2 + [A 2 -(2ZA 7k )2] ' 



(3) 



) ctq (8> cq and the quasiparticle excitation 



where f^2 = a z 

dispersion ^>-k=\J A 2 — (2ZA7 k ) 2 . From the free energy, 
the saddle point equations arc derived 



1 E 

k 



2A 



[2n B (w k ) + l] =3, 



(2^7 k ) 2 
w k 



[2n B (w k ) + 1] = Z. 



(4) 



The dynamic correlation functions of the spin, orbital, 
and spin-orbital tensor operators can be calculated by 
using S? = |*t( r .)nf*( ri ), T<* = i*t( r .)Qa^( ri)) 

* - i*t( r .)0f where Q% = a <g> a <g> a x , 



L 

4 = 



y, iig 



CT0<X>cr C 
- (70<8>cr z 



)o z ; Of, = a C 
? )o~q, and Q^f 



All fifteen correlation functions associated with Sf, T", 
and L" = 2S"T[ are found to satisfy 



(5) 



independent of the indices a and /3, implying the the 
bosonic MF state is SU(4) symmetric! 



2 



On a 2D square lattice and at T — 0, the conversion 
from the summations over momenta to the integrals will 
be invalid as A — > 8 A. By separating the divergent terms 
at k*= 0, Q from the summations, the superfluid den- 
sity p — 4ns(wk*)/-\/l — (8A/A) 2 is defined and obtain 
p ~ 0.107, A ~ 1.3159, and A ~ 10.5271. The dynamic 
correlation function is evaluated, leading to the imagi- 
nary part of the dynamic susceptibility for lu > 0, 



d 2 k 



1 - 



A 2 - (8A) 2 7k7k+c 

£k£k+q 



x [n B (e k+q ) + n B (e k )+l] S (u - e k -e k+q ) . 

As q — * Q with A — > 8A, Im% s (Q, oj) is divergent. From 
the fluctuation dissipation theorem, the dynamic struc- 



ture factor is obtained S(Q,ui) w N 2ir5(u), 
yielding to a long-range AF ordered state with the stag- 
gered magnetizations mg w 0.07566. 

These results disagree with the prediction of a spin liq- 
uid state with a finite gap by the quantum Monte Carlo 
calculations on a 12 x 12 square lattice 0, and is be- 
lieved due to their small system sizes Q. Interestingly, 
the latest quantum Monte Carlo simulations on much 
larger lattice sizes have showed that the convergence 
of staggered magnetization to a finite value can be seen 
for systems larger than 32 x 32 lattices, and the estimated 
staggered magnetization is mg ~ 0.0703, which is nearly 
the same value as ours. 

The author would like to thank S. Q. Shen for his useful 
discussion. 
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